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When Lorentz invariance is violated at high energy, the laws of black hole thermodynamics are 
apparently no longer satisfied. To shed light on this observation, we study dispersive fields in 
de Sitter space. We show that the Bunch-Davies vacuum state restricted to the static patch is no 
longer thermal, and that the Tolman law is violated. However we also show that, for free fields at 
least, this vacuum is the only stationary stable state, as if it were in equilibrium. We then present a 
precise correspondence between dispersive effects found in de Sitter and in black hole metrics. This 
indicates that the consequences of dispersion on thermodynamical laws could also be similar. 
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I. INTRODUCTION 

The two main predictions of quantum field theory in curved space, namely black hole radiation [1-3] and primordial 
spectra in inflation [4-6] share many properties. In particular, both spectra stem from vacuum fluctuations with 
extremely short wave lengths [7]. They are therefore in principle sensitive to the ultra-high frequency behavior of 
the theory. To check this sensitivity, following [8], non-linear dispersion relations, which break the local Lorentz 
invariance, have been used in the context of black holes [9-13] and in cosmology [14-17]. However so far, these studies 
have been conducted separately and with different means. In homogeneous cosmology, time dependent modes with a 
fixed comoving wave vector have been used, whereas for black holes, the analysis was based on stationary modes. In 
spite of this, the two cases are unexpectedly similar, as we shall show. 

In the present work, we analyze dispersive fields in de Sitter space for two reasons. Firstly, since de Sitter endowed 
with a cosmological preferred frame is both homogeneous and stationary, high frequency dispersion can be studied 
along both approaches. This will allow us to relate them in a very precise way. We shall see that their compatibility 
relies on a two dimensional symmetry group which is a sub-group of the de Sitter isometry group [18]. Because the 
generators of the two symmetries do not commute, in each approach only one symmetry is manifest, while the other 
is somehow hidden. In fact, this extra symmetry has been exploited in the black hole near horizon approximation of 
Refs. [9-13], but without noticing (in general) that it relies on properties that are exact in de Sitter space. 

Secondly, the main consequence of high frequency dispersion, that is the loss of the thermality of the spectrum, has 
raised deep questions concerning the relationships between Lorentz symmetry and black hole thermodynamics [19-21]. 
It has even been claimed that this loss should lead to violations of the second law [22-24]. To consider these issues 
in a simpler context, we study the Bunch-Davies (BD) vacuum of dispersive fields propagating in de Sitter space. 
In practice, we analyze the two-point function evaluated in this state. For all dispersion relations, we show that it 
is stationary and periodic in imaginary time with period 2tt/H, where H is the Hubble factor, as it is for Lorentz 
invariant theories. However, in spite of this, we demonstrate that the BD vacuum is no longer a thermal state when 
restricted to the static patch (because the two-point function looses an analytical property which is part of the KMS 
conditions [3, 25]). We also argue that the violations of thermality find their origin in the fact that the stationary 
Hamiltonian restricted to the static patch is no longer bounded from below, and that this loss is not necessarily related 
to the presence of event horizons. (For subluminal dispersion it can occur in stationary backgrounds without Killing 
horizon.) In spite of these violations, for free fields, we show that the BD vacuum is the only stationary state which is 
regular (in a sense that shall be made precise in the text) and that the other regular states "flow" towards this state. 
We expect that this will remain true for interacting theories, such as Xcj)'^. We then argue that the lessons obtained in 
de Sitter should apply to black holes because there is a precise correspondence between dispersive effects in de Sitter 
and in black hole metrics. This aspect is presented in App. B. 

The paper is organized as follows. In Section II we present the basic properties of high frequency dispersion in de 
Sitter space. In Section HI, we demonstrate that the Bunch Davies vacuum is no longer thermal for any (superluminal) 
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dispersion relation. We also show that it is the only stationary stable state. In Section IV, the departures from 
thermality and the S'-matrix are exactly calculated for a quartic superluminal dispersion relation. We summarize our 
results in the Conclusions. In Appendix A, we discuss the group theoretical properties characterizing high frequency 
dispersion in de Sitter, and in Appendix B we study the correspondence between de Sitter and black holes. 



II. DISPERSIVE FIELDS IN DE SITTER SPACE 



A. Ultra-violet dispersion 

We work in 1+1 dimensions and consider the flat sections of de Sitter space. They can be described by ds^ — 
—dt^ + a^dz^, where a = e^* /H is the scale factor, and t the cosmological time orthogonal to the flat sections. 
We assume that the preferred frame associated with high energy dispersion coincides with the cosmological frame. 
Following [21], we describe it by a unit time-like vector field u, treated as a given background field. When considering 
its dynamics it can be shown that u flows to the cosmological rest frame [26-28], like peculiar velocities in expanding 
universes flow to rest. We also introduce the unit space-like vector s that is orthogonal to u. We use covariant 
expressions because we want to be able to transpose the present description to black hole metrics. In the above 
coordinate system {t, z), one has = (—1, 0), and = (0, 1/a). 

The fields u, s define the preferred frequency and momentum by f2 = u^d^, P = s^(}^. The dispersion relation 
then reads 

^ F^{P^) ^m^ + + f{P^). (1) 

We suppose that / vanishes faster than P^ for P — > 0, so as to recover a relativistic relation for P ^ A, where A gives 
the ultra-violet dispersive scale. As such, Eq. (1) can be viewed as the Hamilton- Jacobi equation for the corresponding 
dispersive particle [9, 11]. Using -I- u^u^, — s^s^, this equation reads 

g^'^d^^Sd^S + 111^ + fiis'^d^Sf) = 0, (2) 

where S(t, z) is the action of the particle. On the other hand, Eq. (1) can also be viewed as the dispersion relation 
governing some field. However there is some ambiguity because of the ordering of the differential operators, and non 
minimal couplings. In this paper, we work with [29, 30] 

[-g^'^D^D, + + f {-s^ D ^s'' D ,)] $ = 0, (3) 

where is the covariant derivative. For other approaches based on condensed matter models, see [31-33]. 

At this point it should be observed that the homogeneous Killing field = dz commutes with our u field. In 
cosmological backgrounds, the comoving momentum k = = aP is thus conserved for all dispersion relations /. 
What is peculiar about de Sitter space is that the settings are also stationary. This can be easily seen when using the 
preferred coordinates {t,X) defined by dt = u^dx^^ and dx = s'^^^t, see Fig. 1. (They are often called Lemaitre or 
Painleve-GuUstrand coordinates [34].) The preferred time coincides with the cosmological time, and the new spatial 




FIG. 1: The Penrose diagram of the 2 dimensional de Sitter space. The triangle surrounded by the thick line characterizes the 
Poincare patch covered by the coordinates t,X of Eq. (4). Dotted lines represent t = Cst, and dashed lines X = Cst. The 
static patch — 1 < HX < 1 is the square of the middle of the Figure. 
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coordinate is related to z by X = a{t) z. In these coordinates, the de Sitter metric reads 

ds^ = -dt^ + {dX -v{X)dt)^, (4) 

where v — HX. The vector field Kt — ~dt\x, with the derivative defined at fixed X, is manifestly a stationary 
Killing field. One also verifies that Kt commutes with u. In the presence of dispersion, this is a necessary condition 
for its eigenvalue uj = —dt\x to be conserved. The relation between the preferred and the constant frequency is 
fl = uj-v{X)P. 

Because of these two symmetries, Eq. (2) and Eq. (3) can be analyzed either in the k, t representation as done in cos- 
mology [14-17], or in the w, X representation as done in black hole physics [9-11]. However one cannot simultaneously 
exploit both symmetries because the two Killing fields do not commute. In fact they obey [18] 

[Kt,K,]^HK,. (5) 

This is the algebra of a sub-group of the 2-dimensional de Sitter group SO{l, 2). It corresponds to the affine group of 
M: the algebra contains the translation operator dz and the dilatation operator Hzdz acting on functions of z. Unlike 
the full de Sitter group, it is compatible with dispersion and/or dissipation. As explained in App. A, this follows from 
the fact that both £7 and P are invariant under its action. 

The connection with black hole physics is easily made since the norm of Kt is iff = — 1 + {HXY . It vanishes on 
the (black) horizons located at HX± = ±1. We call them black because the preferred momentum at fixed oj decreases 
as P ~ e~^*, as found near the horizon in black hole metrics endowed with a freely falling frame [9, 11]. In fact, the 
differences between that situation and the present one only arise from the velocity profile v. In de Sitter, v = HX is 
globally linear; whereas, for black holes, v is linear near the horizon region only in a finite domain. As a result, if the 
stationary Killing field Kt = —dt\x is common to both cases, Kz is only an "approximate" Killing field in black hole 
backgrounds, as explained in App. B. 



B. Mode analysis 

Given the two symmetries, the solutions of Eq. (3) can be decomposed either as 



or as 



* = / ^^e'^^^kW, (6) 



_ /27r 

r^e— (7) 

We have introduced the bold notation k to differentiate the norm of k, fc > 0, from k itself which belongs to (— oo, oo). 
In the fc-representation, Eq. (3) gives the second order equation 

J^a„(oa. + -»V^ + /(^))*.«)^o. (8) 

As in all cosmological spaces, see e.g. [17], the general solution thus lives in a two dimensional space and takes the 
form 

Mt)^A^Mt) + {B-kMt))*- (9) 

In de Sitter, and in de Sitter only, the k and t dependences in Eq. (8) can be combined in a single variable which, 
moreover, turns out to be the preferred momentum P = Hke~^^ . Indeed Eq. (8) can be re- written as 

{H^P^dl^^rr? + P^ + f{P^))x{P)=0. (10) 

This possibility is due to the presence of the "spectator" Killing field Kt- Whereas Kz guarantees that k- modes 
separate, Kt tells us that Eq. (8) is invariant under 

k^ke"^, t^t + T, P^P. (11) 

This implies that (j)k{t) only depends on t only through P. The same is true for the action Sk{t) = S{t, z) — kz, where 
S{t, z) is solution of Eq. (2). 
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On the other hand, in the w-representation, the spatial modes obey the higher order equation 

[-{u + tdxv){io + ivdx) - + m2 + = 0. (12) 

Unhke what is found in the fc-representation, at fixed uj, the dimensionahty of the space of solutions now depends on 
the dispersion relation: it is 2n when the highest power of P in /(P^) is 2n. In spite of this it is possible, and very 
instructive, to relate the solutions of Eq. (12) to those of Eq. (8). To this end, it is useful to consider the Fourier 
transform, 

^e-^P^<^„(X), (13) 

where P designates the wave vector, and P > its norm. In the lu, P-representation, Eq. (12) becomes 

[~{uj - iHPdp){Lo - iHdpP) + + P^ + /(P^)] 0^ = 0. (14) 

As in the /c-representation, this is a second order equation (in dp). Moreover one verifies, that 4>i^ exactly factorizes 
as [9] 

4>UP) = P-'^-'^x{P), (15) 

where x is independent of w. In addition, one also verifies that x obeys Eq. (10). These unusual properties are due 
to the other Killing field K^- In fact because of Eq. (11), in {t,P) representation, irrespectively of Lorentz violating 
term f{P^), the modes trivially depend on t through a delta-function: 6{P — Hke^*). When working in the (w,P) 
representation, this implies both the factorization of Eq. (15), and the fact that x obeys Eq. (10). The extra factor 
of 1/P in Eq. (15) is due the Jacobian dt/dP = -1/HP. 

Since Eq. (14) is second order and singular at P = 0, the dimensionality of the space of solutions of (t)uj(P) is 4, 
because P has both signs. The physical meaning of this 4 dimensional space, and its relation with the 2 dimensional 
one found in the fc-representation, are given below. 



C. Scalar product and BD vacuum 



To complete the comparison between the solutions of Eq. (8) and Eq. (12), we consider the conserved scalar product. 
It is given by [8] 

($1, $2) = ijdi (nj$2 ~ **in2) , (16) 

where 11 — —u^df^^ is the momentum conjugated to $. The integral must be evaluated along u^dx^^ = dt = 0, and 
the line element is dl = dX = a{t)dz. 

In the fc-representation, for $k = e'''^(/)k and <I>k' — ^0k', one has 

($k, $k') = 2^(5(k - k') X [a{t) {(bltMk, ~ <\>^'idt(i>l)\ ■ (17) 

The standard normalization ('I'k, ^vj) — 27r(5(k — k') imposes to work with modes (/)k that have a unit positive current 
with respect to a{t)idt. When considering x of Eq. (10), it is convenient to re-express this condition as 

Hxtdpx* - Hx*idpx - 1. (18) 

That is, the x niode is imposed to be of unit positive Wronskian. However, because Eq. (10) is second order, x is 
not completely fixed by Eq. (18). To identify the in mode which describes particles at early time, one has to impose 
that it behaves as the positive frequency WKB mode at early time [2]. Using Eq. (11) to re-express this condition in 
terms of P, the in mode xbd must obey 

g.rF(p-)i4 

XBD - . (19) 



P 

Then the modes 0k with unit positive norm can all be written as 

1 



k(i) - ^ [Ak xbd(P) + (P-k Xbd(P))1 , (20) 
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where Ak and _B_k satisfy \Ai^\ — |i3_k| = 1, and where the extra factor of a/I/Zc ensures that (<I>k, 'i'k') = 27r5(k— k') 
is found when xbd obeys Eq. (18). The state which is vacuum with respect to xbd for all values of k, i.e. -Bk = for 
all k, is the Bunch-Davies (BD) vacuum [35]. 

To handle the mode identification in the w-representation, it is appropriate to work with the Fourier mode of 
Eq. (15) and to separate solutions with positive and negative values of P. In the WKB approximation, positive norm 
solutions describe right moving (U) particles for P > 0, left moving {V) particles for P < 0, and vice versa for negative 
norm solutions. However, the exact solutions of Eq. (14) mix U and V modes. The general solution should thus be 
decomposed as 



H 



AZ Xbd + {B^^ xbd) 



~P) 



H 



(21) 



where the 4 coefhcients weigh the initial (BD) contributions with positive (negative) norm A (_B), and with U ov V 
content. In fact, the scalar product of two such modes 



-iuj t. 



(22) 



This is exact and can be verified be expressing Eq. (16) in the P-representation [11]. Imposing the positive norm 
condition ($(^,$1^') = 2TrHS{uj — uj') on the mode basis constraints the above parenthesis to be unity. Hence, in de 
Sitter, irrespectively of the dispersion relation of Eq. (1), the complete set of positive norm stationary modes contains 
2 modes </'^,(/'^ for lu € (— cx),cx)). One verifies that 2n — 4 solutions of Eq. (12) are not asymptotically bounded 
in X, and cannot be normalized. These modes should not be used when decomposing the canonical field obeying 
Eq (3) [36] . It is interesting to notice that the completeness of the stationary modes follows from the completeness 
of the Mellin transform [37], in a manner similar that the completeness of the homogeneous mode basis follows from 
that of the Fourier transform, see the end of App. A for more details. With this remark, we have verified that the set 
of the asymptotically bounded solutions of Eq. (12) matches that of the solutions of Eq. (8). 

To conclude this Section, we point out that the iS-matrix in the /c-representation factorizes into 2-mode sectors 
containing particles with opposite wave vectors k, because the space in homogeneous. Instead, in the w-representation, 
the S'-matrix factorizes in different sectors with a; > 0, each of them being a 4-mode sector which contains two U 
modes 0^, (0?^^)*, and two V modes 0^, {(fi^i^)* ■ The 4x4 character of the S'-matrix in this representation results 
from the composition of the cosmological mixing of U and V modes with the stationary mixing of modes of opposite 
frequency, see Sec. IV C for details. 



D. The two Hamiltonians 



In preparation for the analysis of the stability of the BD vacuum, we study the Hamiltonian of our dispersive field. 
We first point out that the fields u and Kt define two different Hamiltonian functions, that we call respectively 
and Ht- Using the conjugated momentum H = —u^d^4) and the Lagrangian density L = — <l?0$/2, where 0$ = is 
Eq. (3), they are respectively given by 



Hu^ j dl{Iidtl^-L), 
Ht = J dlindt\x<l)-L), 



(23) 



where we recall that dt\z — —u^^d^ and dt\x — —K'^d^. Hu thus engenders time translations at fixed z, while Ht 
does it at fixed X. In de Sitter space, and Ht differ because u = Kt — vs ^ Kt since the flow v = HX does not 
vanish. In Minkowski space endowed with a Cartesian u field, which is obtained in the limiting case — 0, the two 
Hamiltonians coincide since u ^ Kt when H ^ 0. This implies that in de Sitter Ht and H^ share the properties that 
the Hamiltonian possessed in Minkowski space. 
On the one hand, the stationary Ht 

Ht^2j dtX {{dt<Pf + m202 + (1 _ v^){dx<Pf + 0/(-5i)0) , (24) 

is conserved for all dispersion relations. However, for both Lorentz invariant theories and dispersive ones with 
/(p2) > 0, it is not positive definite precisely because Kt is space like outside the horizons. (Recall that its norm is 
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—K"^ = \ — .) Notice also that when working in the P representation one easily verifies that the last term in Ht is 
positive definite for / > 0. For dispersive theories with / < 0, such as phonons in Helium*, the density of Ht becomes 
negative where v reaches the critical Landau velocity [7]. In any case, in de Sitter, when using the stationary modes 
of Eq. (21), Ht can be decomposed as dujH^, where 

H^=co{\A^^\' + \A^f-\B^^f-\B^j') . (25) 

Ht is thus manifestly conserved and not positive definite. 

On the other hand, the cosmological can be decomposed as H^ = dki/k where 

H^{t) = ^ (^\dtMt)f+F'{^)\Mtf'^ ■ (26) 

Hu is thus positive definite whenever > 0. (Notice that theories with < arc dynamically unstable even 
in Minkowski space.) However, Ht is not conserved because d\na/dt = H ^ 0. The non conservation of H^{t) 
engenders non-adiabatic transitions [38] which describe pair creation of quanta of opposite k, see Sec. IV A for a 
particular example. When using Eq. (20), the time dependence of Hk can be entirely expressed through P = Hke~^* 
as 

^^^^^ = +^^-'-'' (|i/P9pXBD|^ + F^{P') IxbdI^) + Re I {{HPdpXBBf + F\P^) xl^) } ■ (27) 

We also see that when imposing \A]J^ — |_B_k|^ = 1, the minimization of Hk (more precisely, its integral over one 
period) implies i?_k = 0. This is the classical equivalent of saying that the BD vacuum is the lowest energy state 
with respect to the preferred frame field u. 

To conclude, we clearly see the complementary roles played by the Hubble constant H . In the stationary repre- 
sentation, it is responsible for an energetic instability^ i.e. for a conserved Hamiltonian Ht unbounded from below. 
Instead in the homogeneous representation, H is responsible for the time dependence of the positive definite 
which engenders pair creation, i.e. a vacuum instability. These two properties are valid for all dispersion relations, 
and therefore they also apply to Lorentz invariant theories. This is a reminder that field theories in de Sitter space, and 
in black hole metrics, are threatened by dynamical instabilities, i.e. complex frequency modes [39, 40]. In addition, 
as argued below, violations of thermodynamical laws are also related to an energetic instability. 



III. THE CONSEQUENCES OF LORENTZ VIOLATIONS 



In de Sitter space, when considering Lorentz invariant fields, the Bunch-Davies vacuum possesses many remarkable 
properties. On one hand, it is homogeneous and stationary, and on the other hand, it is an Hadamard state. In 
fact, it is the only stationary Hadamard state [41, 42]. In addition, it can be shown that all other Hadamard states 
fiow towards the BD vacuum. By this we mean that the n-point functions evaluated in these states flow towards 
the corresponding one evaluated in the BD vacuum. In this sense, the BD vacuum is the only stable regular state. 
Finally, when evaluated in the static patch \HX\ < 1, the n-point functions are all thermal [25]. They indeed obey 
the double KMS condition; they are periodic in imaginary time with period 27r/i7, and they are analytic in the strip 
< Imi < 2tt/H. 

When considering dispersive fields, we shall see that, for free fields at least, the BD vacuum still satisfies all these 
properties, save the very last. In fact, even though the periodicity in Im t is still exactly found, the analyticity in the 
strip is always lost when there is high frequency dispersion. This means that the BD vacuum is no longer a thermal 
state. 



A. Stationarity and periodicity 

Since we work with free fields and since the BD vacuum is a Gaussian state, we only need to consider the 2-point 
function. When using the settings of the former Section, the Wightman function in the BD vacuum can be written 
as [2] 

/oo ik{z~zo) 
dk XBoikHe-"') x*BD{kHe-"''>) . (28) 
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When considered at fixed to and t, this function is manifestly homogeneous. It is also stationary, when considered at 
fixed Xq = a{to) zq and X = a{t) z. Indeed in terms of P = k/a{t), one gets 

GMX, t, Xo,to) = dP '^^^^^^ '^p''"^'~'°^\ Br>{P) X^D(^e^(*-«°)), (29) 

which is a function of t — tg only. Hence, for all dispersion relations imposed in the cosmological frame, the BD 
vacuum is both stationary and homogeneous, as it is for rclativistic fields. More surprisingly, Gbd is also periodic 
in the imaginary time lapse, with the usual period 2n/H, exactly as for Lorentz invariant fields, and for thermal 
functions. 

We expect that homogeneity, stationarity and the above periodicity will be exactly preserved when considering the 
n-point functions of interacting fields evaluated in the BD vacuum, because these properties are protected by the 
affine group of Eq. (5). In other words, the n-points functions will always be invariant under this sub-group. 



B. Thermality 

For Lorentz invariant theories, it has been shown by Gibbons and Hawking [43] that freely falling observers immersed 
in the BD vacuum detect a thermal bath with a temperature Tff = ^ in natural units fcs=fi- = c=l. It was 
also understood that, when restricted to the static patch — 1 < HX < 1, the reduced density matrix of any QFT 
(interacting or not) is a thermal state at that temperature. Interestingly, this result is always violated for dispersive 
fields. 

To demonstrate this, we consider particle detectors which follow orbits that are stationary with respect to the 
Killing field Kt- Because Kf = — 1 + H^X^, the only stationary orbits are at fixed X, and with \HX\ < 1 when they 
are time-like. The detector transition rate of spontaneous excitation R_ (de-excitation i?_|_) is proportional to [3, 44] 



R±{io,X) = 



Hdt 



GBuiX,t;X,0). 



(30) 



To relate the detector energy gap AE > to the Killing frequency w, one must take into account the X dependent 
red-shift factor {AE = uj/y/1 - H^X^) coming from the detector's kinematics, which also enters in Tolman law 
Tioc{X) = Tgi/Vl — H^X'^ relating the local temperature to the globally defined one [45, 46]. In what follows, we 
shall work with the globally defined temperature and with to. To study the deviations from thermality, it is convenient 
to use the temperature function Tgi(uj,X) defined by 



R^{io,X) 



^-uj/T^iiu.X) 



(31) 



For relativistic fields, one has Tgi{uj,X) = Th = H/2Tr, for all \HX\ < 1 and for < w < oo, in accord with Tolman 
law and the Planck spectrum. To compute Tgi in the presence of dispersion, we shall use the fact that the rates are 
given by 



R±iu;,X)^ (bTj'^iX) 



1,BD,C// 



(32) 



where (t)^^'^ is the positive norm BD mode that is initially right moving, i.e. — 1, = = B^^ = in Eq. (21). 
Similarly, {(i>_J )* is given by the negative norm, negative frequency, mode: B^ = 1. In the above equation we have 
used the symmetry X — > —X to express the contribution of the left moving V-vaode evaluated at X as that of the 
right [/-moving one at —X. Explicitly, (j)^'^ is 



dP 



iPX 



2ttH 



XBU{P), 



(33) 



where xbd obeys Eq. (10). To prove that thermality is violated it is sufficient to work with X — and to consider 
very high frequencies uj/A^l. In this limit the integral is dominated by high values of P, and therefore by the 
leading term of the dispersion relation, that we parametrize here by 



p2r 



(34) 
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In the high P regime, the WKB expression of Eq. (19) offers a reHable approximation of xbd- Hence, up to irrelevant 
constants, one gets 

^2°'^(X = 0) « I ^P-^^-^ e^P\ (35) 
Using Q — P" as integration variable, one obtains a F function, namely 

^ 0) « e-/^"« X r (-^ ^) . (36) 

Using this result in Eq. (32), Eq. (31) gives 

Tgi(w»A,X = 0) = n^, (37) 

i.e. n times the standard temperature Th- 

Instead, for ui/A <C 1 and H/A <C 1, Tg\{uj,X — 0) reduces approximatively to the standard temperature Th [13]. 
Hence the BD vacuum is no longer thermal. This result is non trivial since Gbd of Eq. (29) is still periodic in 
imaginary time, with the standard period. From the above equations and from P ^ e^^*, one understands that the 
power n of Eq. (34) reduces the domain of analyticity of Gbd in Im (t) by a factor of n. Indeed since Xbd ^ e*^ for 
large P, the integral in Eq. (29) contributes as 1/(1 — e"^*) which is analytic in the reduced strip < Im(t) < 2tt /nH 
only. The observation that Tgi(w, X = 0) = nTu for w/A ^ 1 shall be verified, for n = 2, in a exactly solvable model 
in Sec. IV B. 



C. Regular states and stability 

In this Section, we show that some of the ingredients of event horizon thermodynamics [43, 47], are still present when 
adding high frequency dispersion. Namely, we show, firstly that the BD vacuum is the only stationary state which 
is regular and, secondly, that the other regular states flow towards the BD state. (As explained below, the notion 
of regular states should be understood as the generalization of Hadamard states in the presence of short distance 
dispersion.) Hence for free fields at least, the BD state is the only stable state. To prove these claims we shall use 
concepts that are common to Lorentz invariant and dispersive fields. 

Before proceeding, let us discuss our criterion of stability. We say that the BD vacuum is stable because, at large 
time, observables computed in nearby states converge towards those evaluated in the BD vacuum. Hence, for these 
observables, the perturbed states will be asymptotically indistinguishable from the BD vacuum. This flow is often 
referred to as a cosmic no hair theorem [48-50] as it closely follows the Price's no hair theorem [45]. We adopted this 
criterion because there is no stationary Killing field which is globally time-like in de Sitter. As a consequence, there 
is an energetic instability, see Eq. (25), which means that stability cannot be deduced from a spectrum bounded from 
below. It is worth mentioning that to study the thermalization in interacting quantum field theories, the flow towards 
stationary thermal states is established in [51] by studying some n-point functions. Even though the purity of the 
initial state is preserved by the Haniiltonian evolution, after a while, these functions become indistinguishable from 
thermal ones. In that case as well, the stability of the state is thus inferred from the flow of some observables, rather 
than from the evolution of the state itself. 

Since high frequency dispersion modifies the short distance behavior of the 2-point function, we first need to define 
what we mean by "regular states" because the standard definition of Hadamard state is precisely based on this 
behavior [2]. In homogeneous cosmological spaces, this difficulty can be overcome because one can rephrase the 
standard definition in terms of an adiabatic expansion of the solutions of Eq. (8) at fixed k. Since these new terms 
are common to both Lorentz invariant and dispersive field, one can implement the subtraction procedure to dispersive 
fields. Let us recall the key elements, for more details, see [52]. In de Sitter, because of Eq. (5), the adiabatic expansion 
can be done in terms of a single mode x, solution of Eq. (10), and of unit Wronskian, see Eq. (18). This expansion 
generalizes Eq. (19) and is best expressed as 

where W obeys the non-linear equation 

,,,,_ F\P^) (dlW ?,{dpWf\ 

^ ~ P^ 2 [w 2 )' ^^^^ 
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and where determines the dispersion relation in Eq. (1). 

When working with Lorentz invariant fields in D dimensions, the first 1 + D/2 terms in a iterative solution of 
Eq. (39) should be taken into account when determining the 1 + D/2 quantities that need to be subtracted. This 
guarantees that the renormalized stress tensor evaluated in the BD vacuum is finite in cosmological spaces, and thus 
in de Sitter. This is not a surprise since xbd and Xadiab obey the same condition for P —> oo. Hence their differences 
develop at finite P, and because of the expansion H. When working with dispersive fields, this finiteness is still found 
when F'^ is positive, sufficiently regular, and grows faster that P^ for P — !■ oo. Indeed the higher the power n in the 
leading term of Eq. (39), the more suppressed are the next order terms in the adiabatic expansion [52]. For instance, 
in two dimensions, for ^ p2n j ^n-i ^ with n > 2, the second quantity which is usually subtracted in the stress 
tensor is already finite. It can thus be either subtracted or not. [66] In either case, in the BD vacuum of de Sitter, the 
renormalized values oi p — u'^mT^^ and 11 = s^s^Ty^y are constant in space and time, while the flow J = u^s^T^y 
vanishes. 

We now consider the change of the stress tensor with respect to that of the BD vacuum when working with some 
(possibly mixed) state ^' described by the density matrix p^. For free fields, this change is determined by the difference 
of the 2-point functions (5G.p = — Gbd- This difference can be expressed in terms of the positive norm BD modes 



aBD 



as 



^G* = 2Re 



27r 



.(k.-k'.') ^BD (j^^ j^,) (^BD (^,))* ^ k') <^l? (f) } 



(40) 



where n^(k, k') and c^(k, k') are expectation values of normal ordered products of BD destruction and creation 
operators Ok, aj^., : 



n,p(k, k') = Tr 



/5* a^/ak 



c,p(k, k') = Tr [/Sif CkO- 



(41) 



They respectively encode the power spectrum and the coherence of p at the Gaussian level [53] . 

To establish the stability of the BD vacuum, we first point out that the other stationary states are all singular. The 
reason comes from the fact that the stationarity of Gi^ implies that, irrespectively of c*, k x n.p(k, k') only depends 
on the ratio k/k'. Therefore, the change of the expectation value of Flu of Eq. (23) with respect to the BD vacuum, 
necessarily diverges because the contribution of high k is not suppressed enough, as rii^ (k, k) cx l/Zc. This is true for 
dispersion functions F{P) > e > for P — > oo, and therefore true for Lorentz invariant theories. This generalizes the 
fact [42] that the a-vacua, which are invariant under the full de Sitter group and therefore stationary, are all singular, 
save the BD vacuum. 

In our second step we consider states that describe, at some initial time, local perturbations containing a finite 
number of BD particles: TVtot = / cikn,i,(k,k) < oo. Moreover, to be able to handle all dispersion relations at once, 
we suppose that there exists a cut off wave number fcmax above which the number of particles decreases exponentially, 
i.e. 



n^(k, k) < e 



-bk 



(42) 



with 6 > 0. Then Schwartz inequalities and the hermiticity of p<ir implies the following inequalities generalizing those 
of [53] 



|n*(k',k)|' < |n*(k,k)||n*(k',k')|, 



and 



/ 



c?ki/kiC*(ki,k) 



< 71* (k, k) 



y"dkidk2/ki/k.('^*(ki,k2) + <5(ki-k2)) 



(43) 



(44) 



for all test functions /k € C. 

Using these inequalities one can study the behavior of Eq. (40) at large time. Since there is a momentum cutoff 
^maxj at large time only low momenta P matter. Hence for all dispersion relations of Eq. (1), the dominant term is 
the mass term. One should then distinguish massive fields with m > H/2 (see Eq. (48) for the origin of this condition) 
from massless fields. At this point, one also needs to consider the pair creation amplitudes relating the initial BD 
mode ^k° obeying Eq. (19)) to the out mode (/j™* defined at low momentum. Using the techniques of [38] and the 
fact that Eq. (10) is second order for all F^, we can verify that for both m > H/2 and m = 0, the ak,Pk coefficients 
of Eq. (50) are bounded for dispersion relations with F^ > 0. Using this result, at large time and for massive fields, 
one finds that 



(5*G < e 



-H{t+t')/2 



(45) 
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i.e. 6<s,G decreases exponentially in time. This implies that the changes of density, current and pressure with respect 
to the BD vacuum Sp^ = u^^u'^ST^^^ SJ\i, = u^s'^dT^^ and i5n,i, = s^s^dT^fj^ also flow exponentially to 0. On the other 
hand, when m = 0, at large times the dispersive modes become conformally invariant, i.e. proportional to g**'^"''"''' 
where i] oc e~^' is the conformal time. As a result, both scalar derivatives u^d^j^SG^, and s^df^SG^i flow to as in 
Eq. (45). This implies that 5p^, 6J^ and 511^ also flow exponentially fast to 0. 

In conclusion, we have shown that for all dispersion relations, the mean stress tensor Tr(p^ T^^) computed with an 
arbitrary localized state containing a finite number of BD quanta flows towards that computed in the BD vacuum. This 
follows from the cosmological expansion a ~ e^* which redshifts the momenta P ^ fce~^*, and dilutes the particles. 
In our proof we have used the condition of Eq. (42) because, for all polynomial dispersion relations, it guarantees that 
the change of the stress tensor with respect to its value in the BD vacuum is finite. Less restrictive conditions, and 
therefore larger set of states, can certainly be used once having chosen some class of dispersion relations. One could 
also relax the condition that the perturbation is local. However a detailed study of these extensions goes beyond the 
scope of this paper. 

D. Discussion and lessons for black holes 

For all dispersion relations, we showed that the BD vacuum is the only state which is stationary, regular, and 
stable. For Lorentz invariant fields, this stability goes together with the fact that the BD vacuum is a thermal state 
when restricted to the static patch, \HX\ < 1. Therefore, the flow of nearby states towards the BD vacuum can be 
meaningfully considered as a no hair theorem compatible with the laws of event horizon thermodynamics applied to 
de Sitter. For dispersive fields instead, there is a tension because we have demonstrated that the BD vacuum is no 
longer in thermal equilibrium when probed by static particle detectors. 

Before considering possible consequences [22, 23] of this loss, it is important to identify its origin. It can be traced 
to the loss of the positivity of Hf^*", the stationary Hamiltonian operator of Eq. (24) restricted to the static patch, i.e. 
for \HX\ < 1. To explain this, let us first recall the key properties found with relativistic fields [2, 3, 44]. For these 
fields, there exists a complete set of (positive norm) negative frequency modes in which all modes identically 
vanish in the static patch, see Eq. (59) for the massless case. As a result, only positive frequency modes live in the 
patch, and the spectrum of Hf^^^- is bounded from below. 

For dispersive fields, this property is lost because the corresponding negative frequency modes no longer exactly 
vanish for \HX\ < 1. For instance, for superluminal dispersion, they are decaying in this domain, see e.g. [13] for 
details. When Hf'^'^- it is no longer bounded from below, given a certain energy in the static patch, it is no longer 
meaningful to look for the maximum of the entropy. In other words, one is lacking a necessary condition for the 
Ordinary Second Law of thermodynamics (OSL) to hold. In this we do not agree with the claim "Only the validity of 
the GSL is in question" [23], where GSL refers to the generalized second law involving black holes. The fact that the 
violations of thermality are neither necessarily related to black holes, nor to event horizons, is clear when considering 
sub- luminal dispersion relations. In that case, as pointed out after Eq. (24), when v exceeds the Landau critical 
velocity, Ht is no longer bounded from below. As a result, there exist stationary metrics without Killing horizons 
where there will be pair creation of quanta of opposite frequency. In these cases, nothing will guarantee that the 
ordinary zeroth law will apply, i.e. that stationary spectra will be thermal. Similarly, the OSL will not hold if one 
ignores the negative frequency partners, as done in [22-24]. Notice that the importance of the violations will depend 
on the intensity of the mixing of modes with opposite frequency (which governs the instability of the system) . When 
the UV scale A is much higher than the typical value of spatial gradient dxv, this mixing could be strongly suppressed, 
and therefore the violations of the thermodynamical laws accordingly so (because the system will be long living) . 

Having clarified these aspects, we now point out that in de Sitter the non-thermal behavior manifests itself in the 
w-representation, whereas the stability proof heavily used time dependent effects in the homogeneous A:-representation. 
However, we know that these two descriptions ought to be compatible with each other. The lesson therefore seems 
to be that in the presence of ultra-violet dispersion, the BD vacuum is still an equilibrium state, albeit with unusual 
(non-thermal) properties when considered in the stationary picture. So instead of viewing the violations of thermality 
as an indication that horizon thermodynamics might no longer exist, these laws could still exist, albeit in some 
modified and still unknown guise. 

Since the flow towards the BD vacuum of Eq. (45) is engendered by the cosmological expansion, a ~ e^*, it seems 
a priori difficult to apply the stability proof to black hole backgrounds, because these are stationary. However, the 
effects of dispersion in black hole backgrounds turn out to be essentially the same as in de Sitter when the surface 
gravity k = H <^ A (see App. B). Hence, there exist good reasons to believe that what applies to de Sitter could also 
apply to black holes. 

Let us present here the essential aspects of this correspondence. First, when the preferred frame is freely falling 
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(and not e.g. at rest with respect to the stationary KilHng field) the commutator of u and s obeys Eq. (Bl) which is 
the generaUzation of Eq. (A5). Eq. (Bl) in turn imphes Eq. (B4) which gives, near the horizon, P ~ e~'^* both for 
outgoing and infalHng modes, as if they were propagating in an expanding de Sitter universe. Secondly, at the level 
of the quantum theory, the deviations from the relativistic black hole spectrum are governed by the quantity D of 
Eq. (B3) which controls the spatial extension of the near horizon region which can be mapped on a de Sitter space 
endowed with a preferred cosmological frame. This shows that the black hole-de Sitter correspondence is not only 
qualitative, but quantitatively determines the spectral deviations. On this basis, two alternative physical scenarios 
can be envisaged when dealing with a dispersive theory of gravity, such as Einstein- Aether [54] or Horava gravity [55] . 
Either black holes are dynamically unstable in this theory, and there is no question of thermodynamical laws when 
the life time (inverse growth rate) of the instability is sorter or comparable to time scale of the processes under study. 
Or they are stable, and one can conjecture that the modified black hole thermodynamical laws will be quantitatively 
the same as those applying to de Sitter. 



IV. QUARTIC SUPERLUMINAL DISPERSION 

It is of value to explicitly compute the modifications of the observables which are due to high frequency dispersion. 
In de Sitter there are a priori two types of observables, firstly the pair creation rates which are due to cosmological 
expansion, and secondly the thermal-like response of stationary particle detectors. In Sec. IV C we shall study a third 
type of observables, namely asymptotic pair creation rates in the w-representation, which combines the former two 
phenomena. 

To get analytical expressions, we consider the quartic superluminal dispersion, i.e. / — P'^/A^. In this case, the 
general solution of Eq. (10) is given by 

where C and D are the two integration constants, and where M and W are two Whittaker functions, see Chap. 13 



in [56]. For simplicity, we introduced the adimensional quantities /x=y|^ — i, A = A/i7, p = P/H. Using 

Eq. (19) to characterize the initial large P behavior, the unit Wronskian BD mode, when complex conjugated, is 
given by [17, 52] 




A. Cosmological pair creation rates 



To get the pair creation rates, we need to identify the combination of M and W that corresponds to the final mode 
Xout- As Eq. (19) does not offer a reliable approximation for P — 0, the identification should be done using the 
cosmological time t. Using Eq. (8), one finds that asymptotic positive norm solutions are proportional to e~*^-^* at 
large t. When m < ^, is imaginary and the modes grow or decay at large time [42]. Hence it is not possible to 
define asymptotic out modes. When to > there is no difficulty: when re-expressing e"'''^* in terms of P oc e~^*, 
one gets 

Xout (48) 



Using this behavior, the positive unit Wronskian out mode is found to be 



The in — out Bogoliubov transformation is given by 



Xout = ttfcXBD + /3feXBD 



(50) 
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We put a subscript k to the above (fc-independent) coefficients to remind the reader that all these calculations are 
done in the ^-representation. Using Sec 13.1 in [56], one finds, see App. B.2 in [17], 

l^^(A^'^)l'-^2^(l + ^"^x^^")- (51) 

For A ^ 1, up to exponentially small correction, one recovers the relativistic result, i.e. the first term in the above 
equation. For A <C 1, there is an enhancement of the pair creation probability by a factor equal to e^^^. 

Even though the asymptotic out modes cannot be defined for < m < iJ/2, when m = 0, it is again possible to 
define these modes since for t — > oo, they behave as e"*'^'' where drj — dt/a{t) is the conformal time. It is then possible 
to identify the massless out combination of M and W, and to extract the Bogoliubov coefficients. In this case, the 
norm of f3k is 




For A ^ 1, one gets |/3fe| ~ 1/(64A'*), i.e. a power law decrease, unlike what we found above for the massive case. 
Eq. (52) corrects an error in Eq. (131) of [11] but without altering the conclusions of that Section. 



B. Deviations from thermality 



Following Sec. IIIB, our aim is to exactly compute Tg\{uj,X) of Eq. (31) using Eq. (32). To this end, we need to 
evaluate Eq. (33) for quartic dispersion. Using Eq. (47), we get 

i^^TiX))* = /Je^ j^dpe-^^"-p-^^-i W |, (53) 

Surprisingly, it turns out that this integral can be exactly done, see App. C. Since the final expression is a sum that 
converges as 2~" for large n, one can accurately compute the ratio of Eq. (31) in terms of known hypergeometric 
functions. To study the consequences of quartic dispersion, we plot the temperature function Tgi(a;,Ar) in various 
cases. 

In Fig. 2, we plot Tg\{uj)/TH as a function of uj/H, for various values of A, and evaluated at AT = 0, i.e. for an 
inertial detector. Firstly, when w/A and 1/A are both much smaller than 1, we see that this ratio is very close to 1, as 
expected from former analysis [9-13]. In this robust regime, the detector will perceive a Planck law at the standard 
temperature, up to negligible corrections. Secondly, in the high frequency limit, for lu/A^ 1, in agreement with the 
analysis of Sec. IIIB, the ratio goes to 2 irrespectively of the value of A. This last point is not clear from the Figure 
but can be verified analytically from the expressions of Eq. (C9) and the fact that — )■ 1 when A — )■ 0+. 
Thirdly, we see that there is a sharp transition from the robust relativistic regime to a new regime. An examination 
of Eq. (C9) confirms that the transition occurs at a critical frequency Wcrit = A/2. 

In Fig. 3, we plot logj^Q \Tgi(uj) /Th — 1| to study the small deviations from the relativistic regime for u < ujci-a- We 
first notice that the sharp peaks are due to the fact that Tg\{uj) /Th — 1 crosses while decreasing for cj — > 0. A careful 
examination of the envelope reveals that 

|Tgi(w)/rH-l|~e-"^/«+""/4^. (54) 

Hence, at fixed uj, the deviations decrease exponentially with A, whereas, at fixed A, they grow exponentially till lu 
reaches Wcrit- 

In Fig. 4 we study the X dependence of Tgi(a;, X)/Th- This describes violations of the Tolman global equilibrium 
law. We see that the transition from the robust regime to the new regime occurs at different critical frequencies when 
considering detectors following different orbits labeled by X. Interestingly, this dependence can be expressed as 

where ax — H'^\X\/yJl — H'^X'^ is the detector proper acceleration at fixed X. In addition, on the left panel and for 
\HX\ > 0.9, we notice that the low frequency temperature significantly differs from the standard one. This effect is 
related to the broadening of the horizon that was observed in [13, 57]. In those papers, when considering perturbed 
metric profiles v = Wbackgrd + Sv, it was found that the asymptotic black hole temperature differs from the standard 
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FIG. 2: The ratio of T^i{uj,X) of Eq. (31) over the standard 
relativistic temperature Th as a function of ^ , for X — and 
m = 0, and for four values of A, namely 1, 5, 10, and 50. One 
clearly sees that for large values A, the spectrum is accurately 
Planckian and at the standard temperature, until u reaches 
a certain critical value ojcrit, which is equal to HX/2. For 
(jj > ojcrit, T{ljj,X = 0) increases sharply and reaches 2Th- 
This figure is essentially unchanged when we use a massive 
field with < A/2. 




10 20 30 40 



FIG. 3: The logjg of the temperature difference jT(tj)/Tif — 1| 
as a function of uj, for X = and m — 0, and for four values 
of A, namely 10, 20, 40 and 80. We see that \T{uj)/Th - 1| 
increases exponentially in lu until oj reaches Wcrit. It can be 
shown analytically that \T{lu)/Th — 1| follows Eq. (54). 



one when the spatial extension across the horizon of the perturbation Sv is smaller than kx ^ (A/k)^/^. Here we find 
that the temperature seen by a particle detector differs from the standard one precisely when it enters this region. In 
a log-log plot, we have numerically found that the extension of this region (defined by the locus where the relative 
temperature difference is 1%) depends on A with a power equal to 0.675 ± 0.01 in accord with the 2/3 of the above 
references. Two lessons are here obtained. Firstly, the near horizon properties can be probed either by perturbing 
the background metric v, or by introducing a local particle detector, with coherent outcomes. Secondly, since these 
responses are locally determined, they are common to de Sitter and black holes, in accord with the analysis of App. B. 




FIG. 4; The same ratio as in Fig. 2 (on the left) and Fig. 3 (on the right), for m — and A = 50, and for six different positions, 
namely HX = 0, 0.3, 0.5, 0.7, 0.9 and 0.95. On the right, the last two curves have not been plotted since they are too far from 
the other ones. The corresponding values of the acceleration of the detector are a/H = 0, 0.31, 0.57, 1, 2, and 3. One sees that 
T{lu, X) /Th becomes larger than 2 when X 7^ 0. One also sees that the deviations at fixed cj increase with ax- 

Finally it is also interesting to study the behavior of Tgi{uj)/TH when varying A at fixed uj and for X ~ 0. When 
A is large enough, i.e. larger than the critical value Acrit = ^uj/H, the deviations from the standard temperature are 
extremely small, in agreement to what we saw in Fig. 3. Instead, for A — >■ 0, T{u!,X = 0)/Th always flows to 2, with 
a slope that depends on the value of uj/H. An examination of these slopes shows that the slope decrease when lu 
increases: dT/dX\\=Q goes from 1.02 ± 0.005 to 0. This is the behavior at small A. The behavior at large A was given 
by Eq. (54). 



C. Asymptotic S-matrix in the a;-representation 



In this Section, we compute the Bogoliubov transformation between the initial BD modes and the asymptotic out 
modes in the ui representation. In this representation, the modes are identified through their spatial asymptotic 
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behavior, and not their temporal one we used in Sec. IV A. Hence the Bogoliubov transformation can be viewed as an 
^-matrix. This is the description which is appropriate to study the mode mixing on an analogue black hole horizon. 
For more details about mode identification in the w-representation, we refer to [13]. 

In the present case, at fixed w each basis contains 4 modes. Hence the Bogoliubov coefficients form a 4 x 4 matrix. 
This matrix is an element of ?7(2, 2) since the two modes 0^, 0^ have a positive norm, while {(j)^^)* , (</>^(^)* have a 
negative one. In what follows we first study the massless case, and then the massive case m > H/2. In both cases we 
shall see that the S'-matrix possesses unusual factorization properties that are due to the two symmetries governed 
by Kz and Kt . We shall also see that the elements of this matrix combine the cosmological aspects of Sec. IV A and 
the stationary thermal-like aspects of Sec. IV B 

To compute the coefficients of the S'-matrix, we first need to identify the incoming and outgoing modes. At fixed 
w, for quartic dispersion, the general solution of Eq. (14) contains 8 asymptotic branches, 4 for X — ?> oo, and 4 for 
X — >■ — cxD. In addition, when forming wave packets in w, one finds that 4 propagate towards X — 0^ whereas 4 
propagate away from it. The mode identification is based on this second aspect: The 4 incoming modes, are, by 
definition, the 4 solutions that only possess one incoming asymptotic branch. These incoming modes are simply given 
by the Fourier transform of the stationary BD modes (j)^ , thereby showing that the definitions of in modes based on 
their temporal behavior and the spatial one are perfectly consistent. 

To see this, let us consider as an example {(j^^'^)*. Using Eq. (Cll), its asymptotic behavior can be found 
using [58]. Up to an irrelevant overall constant, one finds 

x^±J^ ^ 1) ^ e*(^)-*-t-& + X (Z^)-i+»w-t + (M ^ -Z^)} , (56) 

where x = HX\/\ and where the coefficient Z is 

^ _ 2-'^-'i+'w^(~^/i)^(i-^f +.^) „,,^,^,^^,,,,^,^^ 

^/7F^(i-^f+^|) " ^''^ 

The first term in Eq. (56) describes the incoming high momentum branch, as can be verified by computing its group 
velocity dX/dt = l/d^Pi^, where — dxS^j is the corresponding root of Eq. (2). The last two terms describe the 
4 low momentum outgoing branches. One verifies that they propagate away from the static patch, two for X ^ oo 
and two for X — >■ — oo. In Fig. IV C we schematically represent the space-time pattern associated with a wave-packet 
made with cj)^^'^ . 

We now have to identify the out mode basis, i.e. the four unit norm asymptotic outgoing modes. As in Sec. IV A, 
we treat separately the massless and the massive case. 



1. The massless case 



Since asymptotic outgoing modes have low momentum P, they obey the two-dimensional d'Alembert equation. At 
fixed cj, the equation for the right moving [/-modes is 

{tcj - HXdxW^ = dx4>Z- (58) 
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FIG. 6: In this Figure, unlike in Fig. 1, t = Cst. are horizontal lines, and X = Cst. vertical ones. The two vertical lines 
represent the Killing horizons at HX = ±1. An incoming (BD) high momentum positive norm (7-mode (in thick line) splits into 
four out modes with low momenta: an outgoing positive norm (7-mode (thick line), a negative norm (7-mode (dots), a positive 
norm V-mode (small dashes), and a negative norm V-mode (large dashes). The respective amplitudes of these four outgoing 
modes are given in Eq. (61). To draw these characteristics, we work with m = 0, uj/H = 5 and A/H = 1000. The lapse of time 
spent very close to the horizon is HAt ~ ln(A/a;). It diverges for A — >■ oo, in which case one recovers the relativistic behavior, 
and ultra high initial momenta. 



For > 0, the out [/-modes of positive and negative unit norm are 



lU, out 



= e{i + HX)- 



lU, out \ * 



e{-i-Hx) 



i-l-HX) 



iuj/H 



y/2L0/H 



(59) 



The y-modes 0^'°"*, ((/'^'J"*)* are obtained by replacing X by —X in the above. 

We put the 4 modes in a vector in the following order — {(j)^, {'t'^uj)* , (t>Z ^ both for the BD and the out 

modes, and we define the S-matrix by $5'^ 



5"^^$^"*. We find that 5^ factorizes as 



ak 13k 

al PI 

I3k ak 

K al 







A 

/3S a, 



(60) 



Moreover, the Bogoliubov coefficients ak,f3k are those of Eq (52), and afj , f3^ are the standard relativistic coefficients, 
taken real, and obeying /a^ = e~'^"/^ and |aff p — = 1- To get these real coefficients we chose the (arbitrary) 
phases of the out modes in an appropriate manner. There are only 4 different coefficients in S^j, and they all have a 
clear meaning when considering one BD mode, e.g. 



/BD, U 



lout, U 



iout, V 



B,. 



lOut, V 



y. 



The a^jT^u] coefficients weigh the mode mixing amongst [/-modes of opposite norm, whereas and 
respectively the elastic and the anomalous U — V mode mixing. The norm of these four coefficients obey 



(61) 



describe 



I A, 



\ak\ 



1), 



1/3. 
\B.., 



X (n5 



1), 



(62) 



where = l/^e'^' " — 1) is the Planck spectrum at the standard temperature Th- We see that the deviations 
from the relativistic spectrum are proportional to \ak\ — 1 = |/3fc| ~ A~'^, as those breaking the relativistic U-V 
decoupling. Thus both deviations from the relativistic theory are governed the cosmological pair creation rates at 
fixed k. We also notice that irrespectively of w and A, the elastic \ A^ \ is the smallest coefficient. We finally emphasize 
that these extremely simple results are exact, and follow from the hypergeometric functions 2F2 of Eq. (Cll). 
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2. The massive case 



As in Sec. IV A, the massive out modes should be handled with care. A orthonormal basis for these out modes is 
given by the following right modes (R) 



lOUt 



0{X) 



(Hxy 



lOUt * 



e{x) 



(HX) 



VMHX) 



(63) 



together with the L-modes obtained by replacing X by —X in the above expressions. We have used this R-L separation 
in the place of the U-V one based on the sign of the group velocity, because, for massive modes the asymptotic group 
velocity with respect to the flow v = HX is no longer well defined. 

We now put the 4 out modes in a vector in the following order <i>^ = {^-lJ)* ^ {'t'-uj)*)^ while the 4 in modes 
are ordered in the same order as in the massless case. Defining again the ^-matrix by = 5*^^$°"*, we obtain 



' Oik 








I3k 1 




r 





Ru) 


" 







Pt 





X 





T-^ 










h 












T 





I PI 















R-uj 








(64) 



On the left matrix, the ak,Pk coefficients are those of Eq (51). Hence, as far as this matrix is concerned, we obtain 
the same structure as in Eq (60). Instead on the w-dependent right matrix, the coefficients are 



T 



-iT(ti-uj/H)/2 



•\/2 cosh7r(/_t — (jj/ H) 



R^ 



^1 cosh7r(/i — io/H) 



(65) 



They obey \T^\ + l-R^;! — 1- Hence unlike what was found in Eq (60) the right matrix now describes an elastic 
scattering between modes of the same norm. As a result, the main difference between the massless and the massive 
case is that the final occupation number of massive particle no longer diverge as Th/uj for w — >■ 0. This disappearance 
of the thermal like divergence was already found in [59] in black hole metrics. 



V. CONCLUSIONS 



In this paper we obtain three kinds of results, precise mathematical ones characterizing dispersive fields in de Sitter 
space, those concerning the correspondence between dispersive effects in de Sitter and for black holes, and finally 
more general ones associated with the observation that thermality is violated when Lorentz invariance is broken at 
high energy. 

Concerning the first kind, in Sec. II, we used the group associated with the two residual symmetries of dispersive 
fields in de Sitter to provide precise relationships between the two representations of the field, based respectively on 
the homogeneity and on the stationarity of the settings. The key result is that the homogeneous modes and the 
stationary ones can be all expressed in terms of the single BD mode Xbd(^) a-nd its complex conjugated, where xbd 
obeys Eq. (10) and Eq. (19), see Eq. (20) and Eq. (21). For free fields, all observables are thus encoded in that single 
mode. The algebraic properties associated with the residual group are further explored in App. A and shown to be 
compatible with both dispersive and dissipative effects. 

Having identified this group, we present in App. B the precise correspondence between high frequency dispersion 
in de Sitter and in black hole backgrounds. Because the de Sitter case is also stationary, many aspects are common 
to both cases, with one exception. In de Sitter, when the preferred frame coincides with the cosmological frame, the 
fields u, s obey the affine algebra of Eq. (A5). Instead, in stationary black hole space-times endowed with a freely 
falling frame, u and s obey the local algebra of Eq. (Bl) governed by Q{x), the expansion of u. Since this is basically 
the only difference, the observables of dispersive fields computed in black hole backgrounds, such as the ^-matrix, 
possess the same properties as in de Sitter, up to inverse powers of Z?'^/^A/k, where A is the dispersive frequency, 
K = 00 is the expansion evaluated on the horizon, D = kx gives the extension of the near horizon black hole region 
which can be mapped onto de Sitter, and where the power 3/2 characterizes quartic dispersion [13]. 

In Sec. Ill we show that the two-point function computed in the BD vacuum is still stationary and periodic in 
Imt with period 27r/iJ, as it is for Lorentz invariant fields. In spite of this, we then show that the BD vacuum is no 
longer a thermal state when restricted to the static patch. In particular, we show that the temperature function of 
Eq. (31) is, for ultra high frequency uj/A ^ 1, n times the standard one, where n is the highest power of in the 
dispersion relation of Eq. (1). In Sec. IV B, by considering the response function of particle detectors with different 
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acceleration, we also show that the Tolman law is violated. Even though the BD vacuum is no longer in thermal 
equilibrium, we prove that (for free fields at least) it is still the only stationary, regular, and stable state, as it is in 
relativistic theories [48-50]. In other words, for dispersive fields, there is no (regular) KMS state on de Sitter space. 
We believe that these properties will remain true when considering interacting fields. Finally we explain the origin 
of the violations of thermality in terms of the loss of the positivity of the stationary Hamiltonian restricted to the 
static patch. Whereas this operator possesses a spectrum bounded from below for Lorentz invariant theories, this is 
no longer true for dispersive fields. As a result the ordinary second law of thermodynamics is no longer protected, 
violations of this law are possible, and the system might develop dynamical instabilities. 

In this respect the fact that the BD vacuum is shown to be stable in de Sitter becomes a non-trivial result. Moreover, 
because of the precise correspondence between dispersive effects in de Sitter and in black hole metrics, we conjecture 
that the properties found in de Sitter should apply to black holes, when these are dynamically stable, or sufficiently 
long living, in the adopted theory of gravity that violates Lorentz invariance at high energy. 
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Appendix A: Afflne sub-group and ultra-violet dispersion 

In this Appendix we explore the relationships between the residual affine group based on the two generators Kf , of 
Eq. (5), and the possibility of considering modifications of the field equation that encode dispersive and/or dissipative 
effects. In a sense we are performing a group theoretical approach to dispersion and dissipation on de Sitter space. 
For a similar approach based on the full de Sitter group, we refer to [42]. 

The three generators of 5*0(1, 2) can be taken to be 

Kt = dt-Hzd,_, K,=K_=d,, K+ = -Hzdt+i = j d,. (Al) 

They are linked to the usual generators by K± = ± Ky , and they obey 

[Kt,K,]^HK,, [Kt,K+] = -HK+, [K+,K,]^HKt. (A2) 

We now wish to characterize the set of local differential operators that commute with both Kt and K^- The most 
general local operator acting on scalar fields can be written as 

oo 

0= Yl an,m{t,z)d^dr. (A3) 

n,m— 

Imposing that O commutes with i^^, imphes that the a's depend only on t. Imposing that it also commutes with Kt 
implies an.m{i) — ctn^m e~"^*, where an.m ^ire constants. Hence, O is necessarily of the form 

oo 

0^ J2 (A4) 

n,m— 

where P — —ie^^^dz = —is^df^ is the preferred momentum operator, and Ct — idt = —iu^^d^ is the preferred frequency 
entering in Eq. (1). What we learned here is that the only vector fields that commute with Kt and are the u and 
s fields associated with the cosmological frame. In addition we notice that u and s obey 

[u, s] = Hs, (A5) 

which is the affine algebra of Eq. (5). Hence this algebra is intrinsic to the cosmological frame on de Sitter space. [67] 
Notice also that Eq. (A5) is exact only in de Sitter. However, an interesting generalization of this equation exists in 
black hole spaces, see Eq. (Bl). 
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Our program is to treat O as defining the field equation that generalizes Eq. (3). To this end, we impose that O be 
second order in Cl. This leads to 

6 = ~(f +g{P)VL + h{P). (A6) 

Next we impose the invariance under the discrete parity symmetry 2; — > — z. This implies that g and h are even 
functions of P. In higher dimension, this condition would follow from the requirement of isotropy. The last important 
condition is that O be compatible with a unitary evolution [60]. The proper way to specify this condition is the 
following: the part of O that is even in describes dispersive effects and should be self-adjoint, whereas the odd part 
describes dissipative effects and should be anti-self-adjoint, where the adjoint is defined by 

j Sx^f^^*{d-^) = j d^xy/^{d^<py ^. (A7) 

To sort out the contributions which are due to the expansion, is useful to introduce the self-adjoint operators flga = 
^ (Cl + = Cl + and Cl"^^ = tl"^ + mCl — Then, the "unitary" operators are given by 

= -Clj^ - i (-fsa^se. + f^saTsa) + Fsa, (A8) 

where 7so and Fsa are both real functions of P^. To be more explicit, when applied to e^^^(f>k{t), the field equation 
6$ = gives 

(- dt adt + 2 {isaaY'^-dtiisaaf^ + (F,, + (t>k{t) = 0, (A9) 

where a is the scale factor, and where the argument of jsa and is k'^/a{t)'^. Because of the affine group, this 
equation can be simplified using the function xiP) = (^kit) introduced in Sec II B. In the present more general case, 
one still gets a single equation valid for all (f>i^ modes: 

(n'P'dl - C-^)y^HP^dpC-yf" + {Fsa + ^))x^ 0. (AlO) 

The function F — Fsa + ^ describes the dispersion effects compatible with the affine group, exactly as in Eq. (8). 
Instead, the function ^sa, which multiplies the odd term in f2, describes the dissipative effects compatible it. It 
precisely matches the set of 7 functions introduced in [60, 61] to describe dissipative effects that are local in time, and 
that obey the Generalized Equivalence Principle, which states that the action must be a sum of scalars under general 
coordinate transformations which reproduce those one had in Minkowski space-time endowed with a homogeneous 
static u field. This agreement is non-trivial and follows from the fact that, on one side, the GEP implies that the field 
equation can only depend on the two scalars fl and P defined by the metric g and the u field, whereas, on the other 
side, 57 and P are the two invariants under the generators Kt and of the affine group. 

It is interesting to impose the invariance under the third generator, namely K^. In that case, the only invariant 
operator is O = —O-la + P^ — —Kf + K+K_ + K^K^ which is the Casimir of 5*0(1, 2). We thus see that neither 
dispersion nor dissipation is compatible with the full de Sitter group. We also notice that the affine group has no 
Casimir operator in the sense that the universal enveloping algebra of the affine group has no element, but the identity, 
that commute with the affine group. 

As a final comment, we notice that the affine group is closely related to Fourier and Mellin analysis [37]. When 
working on £^(M), the eigen- modes of —iK^ — ~idz of eigen- value k are the plane waves e*''^, whereas those of 
—iKt — —iH{zdz 4- 1/2) of eigen-value w, are 4>^ — 6{±z){±zy'^^^~-^/^ . The latter live on either side of z — 
and they correspond to Mellin modes. They are complete for cj G M (since invertible) on /^^(M"*"). Hence, to have 
completeness on functions of M, one need two families of Mellin modes, on either side of z = 0, given by 0^. 

Appendix B: Black hole - de Sitter correspondence 

Let us start afresh with a stationary black hole metric and a preferred frame that we describe by a unit time-like 
vector field u. We also introduce the unit space-like vector field s orthogonal to u. To perform the comparison with 
the de Sitter case in meaningful terms we shall use quantities, and coordinates, that are invariantly defined. 
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To get a situation which is closer to that of Sec. II A, we now add two assumptions. We first assume that the 
preferred frame is freely falling, i.e. a'^ = u^D^u^ = 0. This implies that the commutator of u and s obeys 

[w,s] = es, (Bl) 

where <d{x^) = —D^u'^ is the expansion of the u field. This important equation generalizes Eq. (A5). (To obtain 
it, we used g ~ uu + ss which implies —Df^Ui, = u^j^a^ + Bs^Si^, and the Lie derivative Cu{s^ ® s'^) = Cu{g^^) = 
—D^^u" — D^u^^ — uf^a" + u^af^ + 295^5". Hence, [u,s] = 6s + au where a = ^JcFa^.) We then assume that u 
commutes with the stationary Killing field K^. Under these assumptions, when using the preferred frame coordinates 
dr = u^dx^, dx = s^d^, the metric reads 

ds^ = -dr^ + {dx - v{x)dTf, (B2) 

the expansion of m is 8 = dxV^ and the norm of Kt = —dr\x is — — 1 + f^. The location of the Killing horizon, 
where = 0, is taken to be a; = 0. Then the behavior of v in near horizon region (NHR) is v ^ —1-\-k,x, where 
K = 00 is the expansion evaluated on the horizon [62]. When k > and v < 0, one has a black hole horizon, since 
null outgoing geodesies follow x ^ xq e^'^ in the NHR. 

The important lesson here is that under the assumptions of stationarity and freely fallingness, the black hole metric 
and the preferred frame are completely, and invariantly, determined by v{x). Since the de Sitter background fields of 
Sec. HA can be described by the same settings with the extra condition that v^s is linear in x, the comparison of 
dispersive effects associated with Eq. (1) can be easily done for the solutions of both Eq. (2) and Eq. (3). In particular, 
we can already predict that the deviations between de Sitter and the black hole case will be governed by the spatial 
extension of the black hole NHR where v is approximatively linear in x. When v = —1 + D tanh(Kx/D), the extension 
is, roughly speaking, given by \kx\ — D. [68] Using this velocity profile, near the horizon, Eq. (Bl) is given by 

[u,s]^e{x)s^KS (^1-^^+OiKxf^, (B3) 

which clearly shows that the deviations with respect to Eq. (A5) are governed by kx/D. 

We now wish to show that this correspondence is not limited to the background fields, but extends to the dynamics 
of dispersive fields. At the classical level, this is most clearly seen by considering Hamilton's equations. In particular, 
irrespectively of the choice of / in Eq. (1), the time derivative of the momentum p — dxS — s'^c'^S' obeys 

where Xui{p) is the root of Eq. (1) at fixed oj, i.e. with expressed as fl = uj — v{X)P. We learn here that Eq. (B4) 
is the dynamical equivalent of Eq. (Bl). This establishes how the preferred frame algebra imprints the particle's 
dynamics. Having understood that, as long as kx <C D, Eq. (B4) and Eq. (B3) guarantee that p obeys 

p{T)=poe-^\ (B5) 

as in the de Sitter cosmology where P = k/a{t). It is worth pointing out that this exponential redshift applies for 
both signs of p, i.e. for both right and left moving solutions. This correspondence in p-space also applies to the 
classical trajectories in x-space. At fixed w, x(r) obeys w — v{x)p = ±F{p), where p(r) is the solution Eq. (B4), 
and where + (— ) describes right moving trajectories. As long as u ~' —1 + kx furnishes a good description of v, the 
dispersive trajectories x^jir) in the black hole metric are indistinguishable from those in de Sitter, i.e. Xi^{t) is the 
same function as X^^{t) — 1/7? for = k and t = t. 

The correspondence further extends at the level of the dispersive field because the stationary modes still (exactly) 
obey Eq. (12) in the black hole case. Therefore, near the Killing horizon, the black hole Fourier modes (t>uj{p) factorize 
as in Eq. (15), where x "^iH obey Eq. (10) with H = k. At this point we make two observations. First, Eq. (10) 
resulted in de Sitter from the co-existence of Kt and K^, and their algebra of Eq. (5). Secondly, Eq. (10) was used 
in all analytical treatments of the scattering of dispersive modes on a black hole horizon [9-13]. These observations 
raise several questions: 

• What is the relevance of this correspondence for the S matrix ? 

• What is the validity domain of this correspondence in terms of time lapses ? 



• Can we define a field which is approximatively Killing near the horizon ? 



20 



The first question is certainly the most important one. As shown in [13, 57, 63], in the black hole case, when 
A/k ^ 1, the leading deviations from the Planck spectrum at the standard Hawking temperature are governed by 
inverse powers of the parameter D which enters in Eq. (B3). This means that these deviations are in fact defined with 
respect to the corresponding dispersive spectrum evaluated in de Sitter space. This is perfectly coherent because in 
de Sitter, the deviations due to dispersion with respect to the relativistic spectrum are very small, see Sec. IV A and 
Sec. IV C, much smaller than those of the black hole case. In brief, this explains why the parameter D of Eq. (B3), 
which governs the extension of the black hole near horizon region which can be mapped in de Sitter, also governs the 
spectral deviations of the black hole flux. 

Concerning the second question, as far as space is concerned, the validity range of the linearized expression of v 
around = is trivially fixed by D. What is less trivial concerns the lapse of time during which this linearized 
expression can be used, given the dispersion relation of Eq. (1). It is at this level that the separation between the 
background scale k = Qq and the dispersive scale A enters. When A/k ^ 1, the lapse of time during which the 
right moving f7-particles of frequency w ~ k stay in the NHR scales, for quartic dispersion, as kAt ~ \og{D^^^A/K). 
Correspondingly, the accumulated redshift from the high initial momentum till the final one scales as Pin/Pout ~ 
e"^^^ ~ D^/'^A/k. We see that it combines in a non-trivial manner the scale separation and the spatial extension of 
the NHR. In [63] it was explicitly shown that kAt, the adimensional lapse of time spent in the de Sitter like region, 
governs the properties of the black hole spectrum. 

Having clarified these issues, it is worth returning to geometrical aspects by investigating how a vector field = dz 
can be introduced in black hole space-times and to what extent it could be considered as an "approximate Killing 
field" . It should be first pointed out that, a priori, there exist several ways to introduce a new coordinate z. Indeed, in 
de Sitter, K^^ obeys several properties that can be used to define the vector field in the black hole case. For instance, 
the commutator [u, K^^] vanishes. Using this property to define z, one gets the construction of [60] where the black 
hole metric reads ds^ = —dr^ + a^dz^, with a = v{x{t, z))/v{z) ^ e"'*'^ in the NHR. The disadvantage of this choice 
is that the lapse of time during which the exponential is found is much shorter than the lapse Ar we above discussed. 
A posteriori, it turns out that a better choice is provided by imposing that Eq. (5) be satisfied: 

[K,,K,]^kK,. (B6) 

This implies that = e'^^dx is the derivative with respect to the new coordinate z = xe^'^*. We then have the 
following commutation relations 

[Kz,u] = ieix)~K)Kz, (B7) 

and 

Q(x) — K 

Df,Kz^+D^Kzi, = (s^u^ + u^s^). (B8) 

Since k = Q{x — 0), we see that the deviations from the Killingness, i.e. the second equation, and from the 
homogeneous de Sitterness, the first equation, are both governed by the gradient of Q in the NHR, and not from 
Qq — K itself. It is thus geometrically meaningful, and dynamically relevant, to say that a stationary black hole metric 
endowed with a freely falling frame possesses, in the NHR, an approximate homogeneous Killing field obeying the 
affine algebra of Eq. (B6). 



Appendix C: Evaluation of Eq. (53) 

To obtain an explicit expression for {(t>^'^)* of Eq. (53), we shall compute a more general function A{z) to be able 
to exploit some analytical property in z. It is given by 

A{z) ^ e'^'t^+i^/^ / dp e-'P"" W (k, izp^) , (CI) 

Jo Z2 

and which is related to the BD mode by 

((/.^°'^(A))* = e— ("+i)/4./Ae-f^A("+i)/2A(z = 1). (C2) 

V 47r 

To simplify notations, we introduced a = —3/2 — iuj/H, ^ = — iA/4, v — ifJ,/2, x = HX^/X and rescaled p — > 
pVA. Making a rotation in complex p plane of angle 7r/4 one gets 

A{z)^ I dpp^e-P'^^ W{^,iy,zp^) , . (C3) 
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The Whittaker is then expressed as a sum W v, zp^) = B{i/) + B{~v) where [56] 



^' sin27riyr(l/2-iy-e)r(l/2 + i^-0 ^ nW{l + 2,^ + n) ' 



e 

ri=0 

Then the amphtude A is expressed as 

-A{z) sin27rt/ z^A+^iz) - z^" A^„{z) 



TT r(i/2-z/-e)r(i/2 + i/-e) 

where 



(C5) 



*J n— 

-Z. ^.r(i + 2. + n) - ^^'^^ e 



(C6) 



and where /3 = a + 1 + 2i/. The last equahty is vahd only inside the radius of convergence of the power series which 
is \z\ < 1/2. We notice that 2; = 1 is not in the radius, this is why we introduced the extra variable z. Expanding the 
oscillating exponential in x as a series, we get: 



f°° ^pp/^+2ng-p£!;e"'/-* ^-p^ /2 ^ ^ (-l)''x'-e'^^/'^ ^(^_;^^fc)/3_^„^^ /3 + 1 + fc + ! 



2n 

i-l ^ 9 

fc=0 

Using this expression, the sum over n can be done and gives 

"^ (-N/2zx)'^r(l/2 + Or((fc + 1 + /3)/2) 



(C7) 



2fi(^+^-g, ^^^^^ ;l + 2z.;2z) 
Using Eq. (15.3.8) of [56] and Eq. (C5), one obtains 

A( ^ ^ V- (- V2» r(l + (fc + a)/2 + ^)r(l + (fc + «)/2 - i.) 

h r(-e + (fc + 3 + a)/2) ^^g^ 

.1 .1 ^^fc + 3 + a 1, 

2Fr(2+.-e,2-.-e;-e+^^;l-2^) 

The above expressions are all valid for |z| < 1/2. However, since both A and the sum are analytic on C, the result is 
still valid at z = 1. 

It is then convenient to express the 2-^1 as 

2F1 a, 6; c; « = E w ^r^w ^ \ , ^ ' CIO 
1 all 0)1 c + nln! 

n— \ / \ / \ / 

to split the sum between odd and even fc, and to notice that each sum over fc gives an hypergeometric function 
where 



2-F2(l + a/2 + J/, 1 + a/2 - j/; e; n - ^ + (3 + a)/2); ixV2) 



g(a,e)=r(l+ - +;/)r(l + - -zy) " — ^ ,7 ^ — '-^ (C12) 

^ ' ' ^ 2 ' ^ 2 ' r(n-e + (3 + a)/2) ^ ' 

We verified that the above expression is a solution of the 4*'' order differential Eq. (12). From the symmetries of 
Eq. (12), four independent solutions are 

(^^S'^^l^))*, ^T.'^iX), (^^S-'^l-^))*, (C13) 



22 



The last two ones give the V modes evaluated at X. These four functions are independent because they are orthogonal 
to each other when using the scalar product of Eq. (16). 

In addition, to validate our long calculation, we compared the final expression of Eq. (Cll) with the original integral 
of Eq. (53) evaluated numerically with Mathematica® . We found a perfect agreement. 
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